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Abstract 

A quantum field theoretical model for the dynamics of the disoriented chiral condensate is presented. 
A unified approach to relate the quantum field theory directly to the formation, decay and signals of the 
DCC and its evolution is taken. We use a background field analysis of the 0(4) sigma model keeping one- 
loop quantum corrections (quadratic order in the fluctuations) . An evolution of the quantum fluctuations 
in an external, expanding metric which simulates the expansion of the plasma, is carried out. We examine, 
in detail, the amplification of the low momentum pion modes with two competing effects, the expansion 
rate of the plasma and the transition rate of the vacuum configuration from a metastable state into a 
stable state. We show the effect of DCC formation on the multiplicity distributions and the Bose-Einstein 
correlations. 



1 Introduction 

The mechanism of spontaneous symmetry breaking pQ occupies an important role in par- 
ticle physics, cosmology j2] and condensed matter physics [3]. The understanding of this 
phenomenon has lead to cross fertilization between these fields and many new ideas have 
emerged [4f. When a symmetry is spontaneously broken by a vacuum state, it is well known 
that thermal effects at equilibrium can restore the broken symmetry Typically, there 
exists a critical temperature T c at which the effective potential as a function of some order 
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parameter develops an absolute minimum at = and the system undergoes a phase 
transition and relaxes to this minimum at higher temperature. A manifestation of this, 
in a second order phase transition is that at T = T c the order parameter fluctuates at all 
scales and long wavelength oscillations occur 0,13 • For non-equilibrium situations, the dy- 
namics of the transition of a system from a symmetric to a broken symmetry state leads 
to the formation of new structures and the generation of entropy in the form of particle 
productionjHliniEl- in high energy physics, there are two well known examples of symmetry 
breaking: electroweak symmetry-breaking and chiral symmetry-breaking in strong interac- 
tions. One of the main features of QCD, the underlying theory of strong interactions, is 
the spontaneous breaking of its approximate SU(2)l x SU(2)r chiral symmetry. Sponta- 
neous breaking of this approximate symmetry explains the very small pion masses. If the 
symmetry breaking were exact, pions would be massless Nambu-Goldstone bosons. Another 
consequence of this symmetry breaking is the presence of a non-vanishing quark condensate 
in the vacuum. It is believed that at very high temperature a quark-gluon plasma is formed 
in which chiral symmetry is restored, and much effort is being made to explore such a phase 
transition by means of high-energy hadron or heavy- ion collisions 10]. The question is : 
How do we test the spontaneous symmetry breaking mechanism in a non-equilibrium pro- 
cess directly? Can we create a suitable condition so that the vacuum state is disturbed for 
a small region of space-time and observe different excitations and domain structures in the 
vacuum? 

Recently it has been conjectured that this maybe possible for the chiral symmetry- 
breaking in strong interactions given the present high energy physics facilities [TT] EH HSj • 

To describe aspects of QCD related to this symmetry, it is convenient to introduce an 
effective theory, to use a low energy effective sigma model as a model for the QCD phase 
transition because it respects the SU(2) xSU(2) chiral symmetry of QCD with two light 
flavours of quark and it contains a scalar field (E) that has the same chiral properties as the 
quark condensate. The E field can thus be used to represent the order parameter for the 
chiral phase transition. The chiral four- vector of fields (S, 7r), where E represents the quark 
condensate and the three 7Tj are the pion fields. In the physical vacuum, (E, 7r) points in 
the E direction. If the chiral symmetry were exact then there would be a "chiral circle" of 
states degenerate with this vacuum state. In practice, the symmetry is explicitly broken by 
the current quark masses and so there is a unique vacuum. 

Because of this circle of nearly degenerate field configurations, as the chirally restored 
plasma cools and returns to the normal phase, the system could form regions in which the 
chiral fields are misaligned, that is, chirally rotated from their usual orientation along the 
E direction. There has been much recent interest in this phenomenon, which is known 
as a disoriented chiral condensate (DCC). If such a state were formed, it would lead to 
anomalously large event-by-event fluctuations in the ratio of charged to neutral pions. 

A region of DCC can be thought of as a cluster of pions of near identical momentum 
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around zero( coherently produced) with an anomalously large amount of fluctuation of the 
neutral fraction. In order to produce such a state in a quark gluon plasma, the hot plasma 
must evolve far from equilibrium and in particular it must reach an unstable configuration 
such that the long-wavelength pion modes are amplified exponentially when the system 
relaxes to the stable vacuum state. Thus questions of whether a DCC forms and how it 
evolves cannot be addressed in the framework of equilibrium thermodynamics. Techniques 
for applying QCD directly to such situations do not exist at present. To explain these non- 
equilibrium phenomena, we need to restructure the theory of phase transitions to incorporate 
the micro structures (or states) instead of macro structures. For this we need to set up a 
theory at the quantum level. Some attempts have been made [T7] , [TH] , [20] , [2J] > [22] , but to 
our knowledge these are incomplete as they do not incorporate the disorientation of the 
condensate and the evolution of the plasma in one unified picture. In his talk at Trento, 
(23] Bjorken pointed out some of the unresolved issues concerning the DCC called the DCC 
trouble list. They were: 

• (a) Are coherent states the right quantum DCC description, or should one go beyond 
to squeezed states, etc. 

• (b) How does one link DCC thinking to quantum effects in the data, especially the 
BoseEinstein correlations. Is DCC just another way of talking about the same thing? 

• (c) How does one enforce quantum number conservation, especially charge? [25 

• (d) DCC production may imply anomalous bremsstrahlung, due to the large quantum 
fluctuations in charge. Can this be calculated from first principles? 

• (e) Can one really set up the problem at the quantum field theory level? 

In this paper we lay the foundation for a unified quantum picture of the DCC with and 
without orientation and give our answers to these questions. 



2 The Model 



In most treatments of DCC formation, the classic Gellman-Levy lagrangian [29J is used. For 
a more general treatment which encompasses this Lagrangian, we use the field theoretic 0(4) 
symmetric sigma model with Lagrangian density 
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where $j, i = 1..4 are real scalar fields. 

We note here that the potential we use is not the one traditionally used by the early 
practitioners of DCC analyses. Thus, for later comparison with the dynamical equations for 
the pions used as starting points for analysing the DCC, we will later relate the potential 
that we have chosen to the traditional potential in the Gell-Mann, Levy model (22] • The 
Lagrangian density in the Gell-Mann, Levy model is given by 



- v 2 ) 2 + HE. 



The action is 



~(<W^/A-^(<A 



— V 



S = I d 4 x 

The basic object of study is then the chiral field 

0(r, t) = S(r, t) + it ■ 7?(r, t 



2\2 



(3) 
(4) 

(5) 



where S(r, t) and 7?(r, t) are the scalar and vector fields respectively of the 0(4) vector 
= (£, 7r)and r are the Pauli matrices. 

To make contact with the physical constants A, v and H are related to the physical 
quantities, f n = 92Me V (pion form factor) and m n = 138MeV/ c 2 by the relations 
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Note that the fermionic part is neglected here, since the focus is on the condensate formed 
in the symmetry broken phase where quark degrees of freedom are already confined. With 



the transformation, S— > S — v,the minimum of potential is at < <p >=< £ + n >- 



-< 

./V 



v > 2 = f 2 and the usual equilibrium vacuum state is an ordered state < E >=< v >- 
and < 7r >= 0. 

Now consider the dynamic evolution of the system in a hot quark gluon plasma. When 
the temperature T > T c the system reaches a state of restored symmetry < <fi 2 >= 0. If 
the subsequent expansion of the plasma is adiabatic, the< (ft > field gradually relaxes to 
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the equilibrium state as the system cools to below T c . This is called the "annealing" or 
adiabatic scenario. It was first pointed out in reference [T2| that if the cooling process is 
very rapid and the system is out of equilibrium, i.e., in the event of a sudden quench from a 
state of restored symmetry to a state of broken symmetry such as that occurring in a rapidly 
cooling expanding plasma, the configuration of the field will lag behind the expansion 
of the plasma and there is a mismatch of the configurations and their evolution. After a 
quench, the high temperature configuration does not have time to react to the sudden change 
of the environment, thus the vacuum expectation value < <fi > would stay what it was at 
high temperature for a while and then relax to its equilibrium value. This results in the 
formation of DCC domains and is known as the Baked Alaska scenario. 

Another non-equilibrium situation that can arise is one where the system can go through 
a meta-stable disordered vacuum < £ >= f n Cos(9) < n >= f n nSin(6) and then relax 
by quantum fluctuations to an equilibrium configuration. Here 9 measures the degree of 
disorientation of the condensate in isospin space. 

In both these non-equilibrium situations, the canonical approach to the full quantum 
evolution of the fields is extremely difficult to carry out explicitly. However, we may study 
the quantum evolution of the mean field and include fluctuations. In calculating the effective 
Hamiltonian we use the 0(4) symmetric linear sigma model, whose Lagrangian density is 
given by ((TJ). It is easy to see that with the identification <3> 2 = (£ 2 +7r 2 — f 2 ), and, S— > S— v, 
the effective potential in the 0(4) linear sigma model is the same as the Gell-Mann Levy 
model. We will use this identification when considering the evolution of the pion field. 

The choice of metric is dictated by the simple spherical geometry of the problem. Assum- 
ing a homogeneous expansion of the plasma, a Robertson- Walker type of metric is chosen 
with an expansion rate a(t). This will allow us to examine the competing rates leading to 
"rolling down rate" of the vacuum and the expansion rate a(t). The line element describing 
the expansion of the plasma bubble is chosen to be: 



ds 2 = dt 2 — a(t) 2 dx 2 , 

where a(t) is the expansion parameter. 

The action of the 0(4) sigma model in this metric becomes 



(8) 
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where $j, i 



1..4 are real scalar fields. 
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We now use a background field analysis to study the quantum effects. Assume $j has a 
background classical component fa which satisfies the classical equations of motion: 

fl^, = °. in) 

Treat the quantum field fa as fluctuation around classical solution: 

$. _ 0. + 0. (12) 

since fa satisfies the classical equations of motion, 

1 - 5 2 S 

S = S[<f>i] + 2& g$.g$ . l*=^j + ' ' ' • ( 13 ) 

We shall restrict our analysis to quadratic fluctuations only. In addition, we shall also 
drop the term S[(f>i] as it is just a constant additive term to the quadratic action. Therefore 
we shall deal with a quadratic fluctuation action given simply by: 

S2 = \^A§¥}^- (14) 

For the particular scalar field action given above, assuming all fields vanish at infinity, 

j§- = -d^rWj) - - a 3 ^-. (15) 

Imposing the classical equations of motion, we find that 

dV 

d^g^d^i) + a 3 m 2 fa + a 3 — = (16) 

CHpi 

The equations of motion in this metric are 

a ■ -2 1 dV 
3-fa + fa - -^V 2 0, + m 2 fa + — = 0. (17) 
a a 1 ofa 

Since we are interested in the dynamics of the fluctuation field, we shall treat the fluctuation 
field in S 2 as a classical field and S 2 itself as the classical action for its dynamics. While the 
quadratic part in the fluctuations is reduced to: 

S 2 = I d 3 xdt^-{fa- - \(Vfa) 2 - m 2 4>- - fa -l^^fa). (18) 
J 2, a 5§i5§j 
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We can define a Lagrangian density for studying the dynamics of the fluctuations, C, as 
follows: 



Carrying out a Legendre transformation, it is easy to write down the Hamiltonian density 

d 2 V 



where 
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Assume that the fluctuation field <3>i decomposes into its constituents as: 

=< $ > -(f) 



(23) 



V s j 

We have started with an 0(4) quartet of scalar fields in order that we can construct the 
dynamics of quenched pions in the formation of a disoriented chiral condensate. The physical 
fields are defined so that 



7T+ = -^=(tti +m 2 ); 



71- 



-L^-.vr,). 



Analogously, we define the classical background fields as: 



1 / 



v- = -^(vx-ivi), 



following the identification: 

/ «i \ 

V 3 = V 

V ° 

It is easy to see that the action takes the following form: 
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S = f d 3 xdta 3 {tt + 7t_ - 4r(V7r + )(V7r_) - {m 2 + (4\)v + v„ + \v 2 + Xa^n+n^ 
J a 2 



1-2 1 /V7_ \2 1/2 , „. , o\J , \_2\_2 , 1^2 1 /V7^\2 



2a 2 



o*o - ^( V7r °) " g K + (2A)t; + t;_ + 3A^ + Act )7To + -E - — (VE) 
(m 2 + (2A)w + w_ + 3Act 2 + Awf )E 2 - A(w!tt 2 + w 2 tt 2 + 2<TO 3 7r E 
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1 



+ 2t>_t>37r 7r + + 2f + t>37r 7r_ + 2t>_CT7r + E + 2t> + o"7r_E)} 
The Hamiltonian for this action can be written as 

H H neutral H charged H mixed-, 



where 



H, 



neutral 



= 1 



d 3 xdta 3 {^ 



+ ^(Vvr 2 ) + i(m 2 )( 7 r ) 2 + m:: 



+ 



(pi 



2a 6 2a 2 
1 



2a 6 



-7T n 



1 



2a« +^(V2 2 ) + 5 K)(2) 2 + ^(«|-^)(2) 2 } 



H. 



charged 



I d 3 xdta 3 (^{(p 0+ p - 



+ ^(V7T_)(V7T + ) + (m 2 )(7T + 7T_) + \ (fi 2 ± - ^ ± ) (7r + 7T_) } 

#rmxed = J d 3 Xd 3 dt{\{ V 2 + H 2 _ + V 2 _H 2 + ) 

+ (2A)(f + f37T_7To + f_f37T + 7r + CTt^EvTo + f + CT7T_E + CTf_E7T + )}. 

Where we have also put 

O 2 / j2 

" ^ = X[v 2 + 2v 2 3 ] 



a" 



CO. 



A[w 2 + 2w+?j_] 
A[*; 2 + 2ct 2 ] 



and 



(27) 
(28) 



(29) 



(30) 



(31) 



(32) 
(33) 



2v + v_ + v\ + ct 2 = v 2 . 

This is the most general Hamiltonian for the pion sigma system in the background field 
formalism. The background field can now be parametrized through three angles: 



/ vCos(p)Sin(9)Sin(a) \ 
vCos(p)Sin(9)Cos(a) 
vSin(p)Sin(9) 
\ vCos{9) J 



(34) 
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In order to consider all the special cases that are possible in a transparent way, we shall 
simplify the parametrization of the possible form for the background field to two angles, 
9 and p by letting a — |: then, v± = -^Cos(p)Sin(9), v 3 = vSin(p)Sin(9), and 
a = vCos{9). 



3 Quantisation 

Using standard Canonical quantisation techniques the mode decomposition of the Hamilto- 
nian is: 

r d 3 k 1 Ld n t t 

H neutral ~ J ^2 V^"* + a k a k) 

uo Vt 2 

+ ^("1 - 1 )(4 a * + a fc°l + a -fe°fe + a -k a l) 

+^|(44 + 44) + S(^f - x )(44 + 44 + 4 fc 4 + 4,4)} (35) 



charged = J J^y { ~^ V>k b k + C k C k ) 



(27T) 

" + cfcc ^ + 6 - fcCfe + ct - fe ^ )} (36) 



/■ d 3 k \a 3 v 2 cos 2 (p)sin 2 (9) t 

^nixed = J J^y 3 { ^ + &*C fc 

+ c fc c_fc + 4cL fc + c fc 6i + b[c k + 4& f - fe ) 
Xa 3 v 2 cos(p)sin(p)sin 2 (9) 



+ 



(bkCt-k + &fc4 + 4 a fc + c fc°-fe + c I a - fc + c fe 4 + 4 a fc + b W- k 
\a 3 v 2 sin(p)sin(9)cos{9) 



(d k a- k + d k a\ + d[a k + 4 a - fc ) 
+ 



\a 3 v 2 cos(p)sin(9)cos(9) 



(M-k + b k d\ + cj.4 + c fc GL fc + cj.4 fc + c fc 4 + b l d k + b l d -k) } (37) 
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where 



fl 6 " fl 6 " a 6 "^+ a 2^ 



^ = (m| + (38) 



and 



cr a 

,2 



^P = ^ + ^ + A(^ + 2^). (39) 



a" a 

^ ± ( fc ) _ k 
a 6 a 2 



+ m 2 + \(v 2 + 2v + v_). (40) 



An important point to note here is that although Q n (k) and u n (k) are momentum depen- 
dent quantities, for ease of notation we will drop the k dependence for further calculations 
in this section and revive it when necessary for the description of the physical processes. 

It is very interesting to note that if either of i>±,i>3 is zero, then we obtain the usual 
expected dynamics for the pions with back to back correlations. But if we allow for either v± 
or v 3 to be non-zero as may be envisaged in the highly non-equilibrium dynamics involving 
the formation of a metastable DCC, then, we have terms which involve mixed interactions 
of the pions and sigma. 

Clearly, there are two interesting cases to be considered here: v± = and v 3 = (equiv- 
alently, 9 = 0) and the second case occurs with the formation of a metastable, misaligned 
vacuum when p — | . 
CASE 1 

For this case we use v± = and v 3 = (or, equivalently, = 0) showing that H reduces 

to 

d 3 k 



H = J JW { ^ aW + akal) 



(2 

^(Hr - 1 )(4afc + a k a\ + a. k a k + a)_ k a) k ) + ^{d{d k + d k d\) 



4a 3v w 2 1 — * 1 — ' ' 2a 3 



4a 3 



- 
2 

E 



1)(44 + 44 + d- k d k + d}_ k d\)} 



+{^(blh + c k c{) + _ i)( 6 t 6fc + CfcC t + 6 _ fcCfc + c \ r fcfe t )} (41) 



This Hamiltonian has an su(l, 1) symmetry and can be diagonalized by a series of Bogol- 
ubov (squeezing) transformations simply given as follows: In the neutral sector, writing 

A k (t,r) = fi{r,t)a k + u{r,t)al k = U- 1 {r,t)a k U{r,t) 

4(t,r) = u(r,t)a. k + fi(r,t)a{ = U~\r,t)a{U(r,t) (42) 
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A similar expansion for diagonalization is done for the sigma field,with operators D k (t,r') 
and Dl(t,r') similar to the definition of A k (t,r) and A\(t,r) with the d's replacing the a's. 
For the charged sector 

C k (t,r) = nc k + vbl k = U- l (r,t)c k U(r,t) 

Cl(t,r) = f icl + vb„ k = U~ 1 (r,t)cl k U(r,t) (43) 

B k (t,r) = ^ k + vbi = U~ 1 (r,t)b k U(r,t) (44) 
Bl(t,r) = nd k + vb k = U-\r,t)b{U(r,t) 



where 

, = Ca.Mr) = ^[(| + £) + !] -S«M^i[(^ + g)-l] (45) 

and r is the squeezing parameter, is easily seen to give the usual result that fi 2 — u 2 = 1 
for a squeezing transformation. The complete unitary matrix accomplishing the squeezing 
transformation may be written down as 

JJ^ r ^ _ e J ■^^ r ( k : t ){( a l a -k- a k a -k)Hd{d i _ k -d k d-kW ^g-j 

It should be noted here that putting together all our results for the neutral and charged 
sectors, the total diagonalized Hamiltonian is written in terms of various creation and anni- 
hilation operators as: 

1 1 1 

^^{^{(4^ + 3) + (<$Ck + B[B k + 1)} + n E (DjD fc + -)}. (47) 

Since the £ field decouples, we drop all terms associated with it whenever it is not essential 
to our arguments allowing us to write the total dynamical Hamiltonian for the pion fields in 
terms of the observed pion creation and annihilation operators (a, a\c,c^ ,b and tf): 

r d 3 k Q 

H o = J 7^3 it W + ^)Wk + b{b k + 1) + (iv{(c k b_ k + 6 fc c_ fc ) + {c\bl k + b\cl k )} 

+ ^{(fS + v 2 )(a{a k + l) + vfi(al k a{ + a k a^ k )} (48) 
cr 

This completes the analysis for the case when 6 = 0. 
CASE 2 
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The second case of interest is when p — | . The Hamiltonian for this case, reduces to the 
following: 



H, 



neutral 



d 3 k 1 r ui n . + + , 



(2tt) 3 2 a 3 



a¥A | 
2u w 



1 + 2Sm 2 (#))(4afc + a k a\ + a_ fc a fc + al fc 4) 



+ 



We 



a 3 v 2 A 



(44 + 44) + ^ (1 + 2Co S 2 (fl))(44 + 44 + d- k d k + dl k dl)} (49) 



2w s 



H charged 

H m ixed 

where (for p — |), 



d 3 k t ! (j n 

7^{^(b{b k + c fe 4) + ^(^ - l)( 6 t 6fc + CfeC t + h _ kCk + c t_ fe6 t )} (50) 



rf 3 i AaV 2gm(g)Cos(g) 



(2k) 



) [d k a^ k + 44 + 4 a fc + 4 a -fc)} (51) 



f] 2 


k 2 


a 6 


a 2 




k 2 


a 6 


a 2 



+ mi + Xv z 



+ ml ± + \v 2 (l + Sin 2 (9)) 



(52) 



Unlike the case considered above, we now have a non-zero mixing term coming from the 
7To-S sector. Since both of these are neutral sectors, we can combine the H neutra i and H mixe( i 
terms to form a single neutral sector Hamiltonian which we will again call H neutra i while 
H c harged remains unchanged. 



H, 



neutral 



d 3 k 



(2tt) 3 

I / a-k 




dk 



4 V aP 
Xv 2 a 3 



4 ^ a 6 

^Sin(29) 

Xv 2 a 3 



k_ 



, Sm(20) 
/ ^(l + 2Sin 2 (6)) 



4 



+ 



4 
4 



Sm(20) 
{l + 2Sin 2 {6)) 



Sin{26) 




153) 



These terms clearly show the mixings between the forward and backward pions as well as 
the mixing between the n and E fields. The misalignment of the vacuum through an angle 
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6 induces a mixing of the two fields. The mixed fields are 

( Mk) \ _ ( Cos{9) Sin(9) \(a k \ 
{ D e(k) J \ -Sin{6) Cos{6) ) \ d k j 

Then H neutra i can be expressed as: 



Here, 
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The diagonalization procedure for this case along with the dynamical consequences of this 
Hamiltonian will be discussed in [ST]. 
We finally get 

Hneutrcd = Tfl\ {A\ {k) A e{k) + Aq^A^^) + 71 1 (Aqq,) Ag(_ k ) + >4(fe)^(-Jfc)) 

+ m 2 (Dl {k) D e{k) + D e[k) D ] m ) + n 2 (p e{k) D e{ _ k) + £)J (fe)J Dj ( _ fe) ) (60) 
We now apply two squeezing transformations 

F k = pA e(k) + vA\ { _ k) 



Fl = pA\ {k) + pA B{ - k) (61) 



and 



G k = pDg( k ) + aDl(_ k y 
t _ „nt 



Then, 
and 



G\ = P Dl {k) +aD e{ - h) . (62) 
[**,*H]=1 (63) 



[G k ,G{] = l (64) 
implies that the transformations are indeed Bogolyubov transformations as 

p 2 - v 2 = 1 (65) 

and 

p 2 -(T 2 = l. (66) 
The squeezing parameters p,v,p,a are given by: 




+ 1. (67) 




2 _ jtt n (tt n 1 /t^ ^ / 1 

" V0 S UA 2\LJ + aA i+ 2Vo; s< 

+ \Z^(^( 1_ 2\/^) + ^( 1 + 2 \/^)) 

- 1. (68) 
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2 n w (n* ( 1 /c^\ w,/ i /^v 

2p = few 1 + 2fc 1_ 2fc. 



+ J^(^fi + -^U^fi- 1 ^ 



f^7r \ ^e V 2 V c^tt / fi^ V 2 V 
+ 1. (69) 




+ V^&( 1 + 5^) + S( 1_ 5^)) 

- I- (70) 



With these definitions the neutral sector Hamiltonian is simply: 

V^TT^E ( ( 771 . 1 \ . /V,f ^ , 1 



#ne^ = (^F fc T F fc + - j + + - ) ) (71) 



Combining with the charged sector 



Hcharged = —^{C\C k + B\B k + 1} (72) 

a 



the total Hamiltonian for the case p — f 



tf*/2 = / (0a {^(^ + 4** + 1) + ((^ + ^) + + I))} (73) 



with 



(7 2 k 2 



a 6 a 2 

^ = ^+m^ ± + A^(l + ^n 2 W) (74) 

This completes our analysis of the quantization of the Hamiltonian for the two cases 
mentioned above. We have taken an 0(4) sigma model and succeeded in quantizing the 
quadratic fluctuations to arrive at two Hamiltonians which provide all the required ingredi- 
ents for analysing the formation, evolution and eventual decay of the DCC. Indeed, in case 
2, the Hamiltonian shows explicit mixing with an angle (0), providing a mixing parameter 
and therefore a misalignment parameter in isospin space. 
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4 Evolution of the fluctuations and parametric amplifications 

In the last section we have constructed the quantum Hamiltonians for two cases (6 = and 
P = § )of the 0(4) sigma model. We have also shown the diagonalization of the Hamiltonians 
so that they can be written in terms of the appropriate quantum fluctuation fields as purely 
quadratic Hamiltonians. In this section we shall explicitly consider the case when 6 = (the 
dynamics of the case p = | will be given in a subsequent paper [3T] ) . We notice that H has 
the form of a decoupled Hamiltonian. This is easy to understand from the >S0(4) parent. 
The SO (4) vector has been decomposed into four fields: 7T±,7To and E being respectively the 
charged pions, the neutral pion and the sigma fields. 

It is the sigma field which decouples in this particular Hamiltonian and therefore it can be 
analyzed independently of the pion fields. As our interest is in the pion fields, we can write 
the total dynamical Hamiltonian for the pion fields in terms of the observed pion creation 
and annihilation operators (a, a' , c, c', b and &t)in terms of the squeezing parameters as: 

r d 3 k £1 

H = J J^p-jt W + " 2 ){cW + b{b k + 1} + pv{(c k b^ k + b k c_ k ) + (4&L fc + b\cl k )} 

+ -t{(p 2 + v 2 ){a} k a k + 1} + (vp){a ] _ k a\ + a k a_ k }} (75) 
cr 

We now define the bilinear operators 

V = a k a^ k + b k c_ k + c k b_ k = K{ + K% + K% 
2?t = al k a\ + cl k b\ + bl k c\ = K+ + K+ + K+ 

N = ~{a\a k + a ] _ k a. k + b[b k + bl k b- k + c\c k + d k c. k + 3} = K\ + K° 2 + K% (76) 

and it is easy to see that they satisfy an su(l, 1) algebra 

[N,V] = -V; [N,T)i] = P f ; [V\V] = -2N. (77) 
The su(l, 1) invariant Hamiltonian for the pion fields assumes the form 

r d z h 1 

H= / — ^2n n (k,t)(p 2 + v 2 )N + 2n„(k,t)^(V + V^ (78) 
J (2tx) 6 a 6 

The time dependent evolution equation for the eigenstates of H is given by 

H(t)\m>=^ t \m> (79) 

The particular su(l, 1) structure elucidated above provides us the solution: 

\^(t) >= J T^ rk{v *- Vk) \^(0 > (80) 
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for the evolution of the wave function immediately jooj. Here r k is the squeezing parameter 
related to the physical variables Q n (k,t) and ^(fc) through 

Tanh{2r k ) = f (81) 

Where ft n (k, t — > oo) = u n (k). Thus in the evolution of the condensate , it is the frequency 
changes which bring about squeezing 

The diagonalized Hamiltonian Hq can be converted into a Hamiltonian in terms of quan- 
tum fields corresponding to the operators A, B, C and their adjoints to obtain a purely 
quadratic Hamiltonian (the starting point of many early works on the subject of the DCC) 18, 

! 

We can, for example, write 

^(A[A k + h = §(4A + A k A\) = t) + Pl A (k, t). (82) 



Similarly for B and C, we write: 

<X (BlB k + h = %B\B k + B k B\) = (%) 2 Hi(M) + t) 



cr I cr a 6 

^(ctc k + h = %c\c k + c k ct) = (§) 2 n^(M) + Pl c (k,t) 

(Jj £i CL (Jj 



The Hamiltonian Hq can then be written as: 

d 3 k ^ l„fi 



(83) 
(84) 



H °®=fw^ E k(^i(k,t)+Pl(k,t)) (85) 

J l Z7r J i=A,B,C A a 

The Schroedinger equation for each momentum mode is simply: 

H (k,t)1>(k,t)=i^l>(k,t). (86) 

If we use the n-representation (co-ordinate space representation) for ip(k,t), then, the 
su(l, 1) symmetry of the Hamiltonian tells us that the solution for ip(k, t) is just a Gaussian. 
For simplicity we work with a Gaussian of the form 

< 7r ,7r + ,7r_|^ > (t) =Y[Li(t)exp^ w *^ (87) 

k.i 



while the complete wavefunction is: 



H 3 k 

r(t)= I (2^(M)- (88) 
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Then for each mode A,B,C,D the wave function ip kt i(t) evolves as 

.9^(0 = ( M _ iU 2 {k)Wi{K f))Mt) m 

While 

Combining the above we find that Wi(t) and ipk,i(t) evolve as: 

W) = "»V2|f (91) 

a ■ n 2 

Mt) + 3-Vi,fc(*) - Hr^W = °- ( 92 ) 
a er 

The equation satisfied by the wavefunctions for each mode are then given by: 

Mk, t) + -ip A + (%) 2 (A;, t)MK t) = 
a a 6 

$ B (k, t) + ——IpB + {^) 2 {k, t)lf> B (k, t) = 

a a 6 

MM) + —4c + (%) 2 (M)lMM) = (93) 



where 

(§) 2 (M) = £)+™l + \v 2 (94) 

Where it may be recalled that A k (t), B k (t) and Ck(t) are related to a k ,b k and (the 
physical pion operators), by the squeezing transformation given by 

A k (t,r) = U-\r,t)a k U(r,t) (95) 
B k (t,r) = U-\r,t)b k U(r,t) (96) 
C k (t,r) = U-\r,t)c k U{r,t) 

Cosh{r) = ^/I[(jk + £) + !] v = Sl nh(r) = ^[(^ + g - 1] and r is the 

squeezing parameter. 

£/( r ^ = e / ^^ r ( fc . i )(( a l a lfe- a fe a -fe)+4 d -fc- d fe d -fc)+( c fe 6 -fc+ b fc c -fc)-(4 6 -fc+ 6 l c -fc)) 

The expectation values of the number operator for the neutral pions for each momentum 
k is given by: 

< Mt)W\a k \Mt) >= Smh 2 (r) =< ^ k \Al(t)A k (t)\^ k > . (98) 
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A similar expression may be obtained for the charged pions. 

Since we are considering the expansion of the plasma in a spherically symmetric manner 
as emphasized by the Robertson- Walker type metric ds 2 = dt 2 — a(t) 2 da?, it is possible to 
scale the time so as to provide a conformally flat metric: we let 

d v = a{ty l dt (99) 

so that 

ds 2 = a{r]) 2 {dr] 2 - dx 2 ). (100) 

The equations of motion given above can be transformed into ones that resemble a harmonic 
oscillator with time dependent frequencies. We shall write only the generic form of the above 
equations: In terms of the scaled time, 77, we have: 

In' 

r + —i>' + (k 2 + (ml + A(< $ 2 > -fl))a 2 ii = (101) 

(X 

where a prime denotes differential wrt rj. Here we have a few remarks. Firstly we note 
that to make contact with the dynamical equations for the pions used as starting points for 
analysing the DCC, we return to the point about the potential that we have chosen vs the 
traditional potential chosen in such studies. The traditional potential is ($ 2 — f 2 ) 2 while we 
have taken a potential ($ 2 ) 2 . Therefore, the relationship between our work and the early 
studies is accomplished through a replacement v 2 — > (< $ 2 > (rj) — f 2 ). This explains the 
equation written above and the further analysis of this paper. Lastly, let us scale ip: 

f = aij (102) 

so that the equation becomes: 

- e + vm = (k 2 + mi )e (103) 

where 

V(V) = a- 1 ^ + m 2 n (l - a 2 ) - A(< $ 2 > -f 2 ). (104) 
Thus in the symmetry broken stage < $ 2 >= f 2 

VM = a- 1 ^ + m 2 (l - a 2 ) (105) 

And in the symmetry restored stage< $ 2 >= 

VM = a- 1 ^ + ml (1 - a 2 ) + Aa 2 / 2 . (106) 



19 



These equations have a dual personality: on the one hand they are Schroedinger like 
equations with rj corresponding to the "spatial" like variable and E = lo 2 . Therefore, they 
allow calculation of the reflection and transmission coefficients over the "potential barrier" 
provided by the V{r\) term. On the other hand, they can also be looked upon as equations 
for time dependent harmonic oscillators with time dependent frequencies given by Q 2 and 
uj 2 . These two pictures enable us to calculate the squeezing parameter dependent number 
operator N(k) = Sinh 2 (r(k)) involved in the evolution of the plasma as can be seen in 
appendix 1. 

We also note that the expansion coefficient a(rj) provides us with a control parameter on 
the expansion rate of the plasma while the system provides us with another parameter, which 
we call r for the "rolling down" of the fields in the potential. According to the pioneering 
papers EJ], when the rate of expansion is greater then the relaxation rate of the field 
the DCC forms. 

We now consider two cases: the first being a toy model [33] without the effects of expansion 
and the second being a more realistic "Baked Alaska Situation". 
In the first case, let us suppose that a{rj) = 1 so that 

V(rj) = -A(< $ 2 > -ft). (107) 

Then the equation for £ is: 

t" + (k 2 + m 2 + A(< $ > 2 -fM = 0. (108) 

If r is the time the state spends in the symmetry restored phase, a quench can be modelled 
by assuming that for — | < 77 < | the vacuum expectation value < $ 2 >= and for 77 > ~ 
the vacuum relaxes to its value < $ 2 >=< v 2 >~ f 2 ( since the potential is translationally 
invariant, this could easily be mapped to the interval {0,r}, but our choice makes the 
potential symmetric) . For this sudden quench the problem reduces to that of transmission 

2 

through a symmetric rectangular potential barrier of height (A/ 2 ) = and width r. 

The transmission coefficient for such a barrier is easily calculable ^2J E] and in terms of 
physically known values is: 

T = ^ " • ( 109 ) 

l + ( 2 { t )2 )Smh 2 [{^-ujl)lT] 

From appendix 1 we get 

1 — T 
N(k) = — = 



Sinh 2 [(^-u 2 k 



U k\^2 W k) 



2 t] 



110) 
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The dependence of N(k) on k for different values of r is shown in fig. 1, clearly exhibiting 
the amplification of the low momentum modes. 

NHkL 




Figure 1: Shows the variation of N(k) and k{ in units of m^) with r = .8, .9, 1.0 in units of m % for the 
quenched limit 

Figure 1 also shows us that the longer the system stays in the state of broken symmetry 
the larger the DCC domains. The dependence of N(k) on r for different values of k is shown 
in fig. 2. Recall that by definition, the amplification of the zero modes constitutes DCC 
formation. Since we have a " Schro dinger" wave equation that is exactly solvable, we can 
also calculate the size of the DCC domains. We shall leave that for a subsequent paper |31j . 

This case is similar to that considered by ref [33J. 

For a more realistic scenario, the expansion must be included to show that the enhance- 
ment of the low energy modes and the squeezing parameter are dependent on the rate of the 
expansion mechanism by which the symmetry is restored. To produce substantial squeezing, 
we require a quenched scenario. To show this we have to compare the situation of a sudden 
quench with a slow adiabatic relaxation of the system from the symmetry restored stage 
to the symmetry broken stage. The transition between quenching and adiabaticity can be 
modeled in two ways in view of the dual nature of eqn (jlU4|) . 

We consider the expansion coefficient a(r]) to be of the form 

a 2 ( V ) = ei-^Tanhi—ir] + T -)) + Q( V )Tanh[^( V - £)]. (Ill) 

Here t measures the rate of the expansion and we choose b to be in units of m n . When 
br » 1 we have the quenched limit (fast expansion) and enhancement of low momentum 
modes should occur (DCC formation). When br << 1 we have the adiabatic limit and no 
enhancement of low momentum modes should occur (no DCC formation). 
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Figure 2: Shows the variation of N(k) with T(units of m„ 1) for k = 0, 1 and 2 (units of m w ) for the quenched 
limit 

When viewed as a time dependent oscillator equation in an expanding metric we may 
write the equation for £ in the symmetry restored phase as an oscillator equation 

£" + a4£-V r (77)£ = (112) 

and in the symmetry broken stage as 

£" + u&-V 6 (77)£ = 0. (113) 

Where V& (77) and ^(77) are given by (|l(J6ll(J5j) with the particular choice of the expansion 
parameter a(r]) given by (Jlllj) . Thus the change in frequency from the restored to the 
broken stage is given by 

V r (v) ~ V b ( V ) = +^a( V ) 2 (114) 

Given that we have to calculate the transmission coefficients we need to look at the wave 
functions £ in the limit rj— > ±00. This is possible by solving the oscillator equation with 
a variable frequency: 

r+^ + ^[i-« 2 (^ = o. (us) 

We see that the above equation ()115|) satisfies the required limits fl(r]) 2 = uj 2 as rj— > ±00 
and between these limits it gradually goes to a maximum at = Q 2 . The equation ()115|) 
can be converted into a Schrodinger wave equation for a Wood -Saxon potential barrier given 
by 
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v(rj) = v (e(-,)( i + e ! 6(;?+i) ) + efox^pi^).) (lie) 

Where Vo = ^p. The comparison of the rectangular potential barrier with this barrier is 
shown in fig.3, which reveals it to be a good approximation in the adiabatic limit. 




Figure 3: Shows the variation of the Wood Saxon Potential Barrier of width r = 4 with b. Large br gives 
the Rectangular Potential Barrier (quenched limit) 

We take the values of the variables E = (k 2 + uj 2 )^ and k' 2 = E — V . The transmission 
coefficient for this barrier is given by: 



Sinh 2 (nf )Sinh 2 (ir^f) 
Stnh 4 (ir^±){4\C\Sin 2 (k'b) + (\C\ 2 - 1)} 



T ™ ~ 77. , „, .(k-k')s 77, o/,„n b /i^io 777 ', E > V o (H7) 



Sinh 2 (7T 2(k+k ' ) ) 

and the same with k > %k for E < V where C = m^wr - 

Sinh 2 (ir -^-g — - ) 

Note that for b large this reduces exactly to the rectangular potential barrier (quenched 
limit) transmission coefficient and for b very small this goes over to the Poschl- Teller (Eckart) 
Barrier transmission coefficient. 
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The number of particles of mode k equals 



w(h , _ 1 - (Sinh\n^){4\C\Sin*(Vb) + (|Cf - 1)}) 
[) ~ Sinh?(*%)Sinh?(*%) [ } 

Here b measures the duration of the quench. N{k)vs.k is plotted in fig. 4 . We see that 
in the adiabatic limit of small br , N(k) is exponentially suppressed so that there is no 
enhancement of low momentum modes. 
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Figure 4: Shows the variation N(k) with k for values of b in the adiabatic limit for Wood- Saxon barrier . 

Figure 5 shows the variation of N(k) with k for large values of b , showing the enhancement 
of low momentum modes in the quenched limit. 

From the above we conclude that if expansion time (r)is faster then the rolling time (r) 
we get the quenched limit, while, if expansion time is slower then the rolling time we get the 
adiabatic limit. Since in the squeezed state description N(k) = Sinh 2 rk ,where = r(k) is 
the squeezing parameter , we conclude that in the quenched limit the squeezing parameter is 
much greater than in the adiabatic limit. This demonstrates clearly the connection between 
the rate of expansion and squeezing and the formation of the DCC. If squeezing is large we 
have formed a DCC , if not then there is no DCC. This will now enable us to give signatures 
of DCC formation which are related to the formation process of the DCC. 

5 Pion Radiation from DCC 

Having shown that the dynamics of the evolution of the DCC suggests a squeezed state 
treatment of emerging pion waves and that this effect is more pronounced in quenched 
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Figure 5: Shows the variation N(k) with k for values of b in the quenched limit for Wood- Saxon barrier . 

scenario than an adiabatic one, we now proceed to show the possible experimental signals 
that would result. A complete treatment, with the incorporation of isospin and disorientation 
in isospin space for arbirary momenta is given in jSU] and j3T] . For the present however, since 
we see an amplification of the low momentum modes we consider only the case when k = 0. 
In such case for pions near zero momentum k — > the state (JHUj) factors into a squeezed 
state for the neutral pions and a Caves- Schumaker state for the charged pions. In this limit 
the bilinear operators of equation (fTUj) are 

V = a 2 + b c + c b (119) 
pt = al 2 + clbl + blcl 
N = a\a + +&J&o + ejeo + 3/2, 

Thus the k=0 wave function is 

|^ >= e^-^lV^O) >= e ro( ^ 2+c i^ + ^ c ^ a o +6oco+co&o) |^ (0) > (120) 
The pion multiplicity distribution is given by 

Pn on+ n. = I < n Q ,n + ,n_\^j > | 2 =< n |e ro(a o) 2 - ro * a o|0 >< n + , n_|e 2ro(fe ^- fe ° co) |0 > 

(121) 

defining S(r Q ) as the one mode squeezing operator 



S(r ) = < n \e roi{a ^ 2 ~ a ^\0 > 
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S tm (r ) is then the two mode squeezing operator 

< n + ,n_|e (ro(b £ c £- 6oco)) |0 > 

= S5U_,o- (123) 

The neutral and charged pion distribution is: 

Pn ,n c =< S n()fi > 2 < 5' t ™ + ,„_,0,0 >2 ( 124 ) 

which is just the product of squeezed distributions for charged and neutral pions and only 
even number of pions emerge. Writing n + = n_ = n c , we get the distribution of charged 
particles to be 

_ (tanh(r )) 2 ^ 

nc "tT °' nC ~ (cosh(r )y [iZb) 

and 

_ v _ n !(tan/.(r ))"Q 

° ~ tT ^ ((? y.) 2 cosh(r )2"° 1 DJ 

The generalized squeezed eigenstate leads to products of two types of squeezed states of 
pions at zero momentum. The neutral pions being in a one mode squeezed state and the 
charged pions being in an SU(1,1) coherent or two-mode squeezed state. Thus the neutral 
and charged pion distributions are significantly different as the two types of states have 
different properties. We now illustrate the effect of squeezing in these two distributions. 
Figs. 6, and 7 show the difference in the charged and neutral pion distributions as we vary 
from the squeezing parameter from a low value to a high value. 

We now illustrate the effect of quenching versus adiabaticity on these two distributions. 
Figs. 8 shows the difference in the charged and neutral pion distribution for as we vary 
from the adiabatic limit where the difference is negligible to the quenched limit where the 
difference is significant. 

We now calculate the correlation function by first calculating the the generating function 
corresponding to the multiplicity distribution P n . It is given by Q(X) = Sa(1 ~~ X) n P n . The 
two particle zero momentum correlation function G 2 (0) is given by 

2 _ <n 2 >-<n> _ @|a=i , 
G " <n 2 > " (fj^ (127) 

the variation of G % eutral with the squeezing parameter is given in figure 10. 
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Figure 6: Shows the variation of P nc with n for the ro = 3,3.5 and 4 
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Figure 7: Shows the variation of P„ c with n for the ro = 3,3.5 and 4 
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Figure 8: Shows the variation of P n o (solid line)and P nc (dashed line) with n for the adiabatic limit (r 
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Figure 9: Shows the variation of P„o (solid line)and P nc (dashed line) with n for the quenched limit (r 
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Figure 10: Shows the difference in correlations of the charged(dashed line) and neutral pions(solid line) as a 
function of the squeezing parameter ro 

Thus we see that in the adiabatic limit G 2 neutrol ~ G 2 charged1 whereas in the quenched limit 
G 'neutral <« G 2 charged1 giving a very clear indication of the effect of DCC formation on the 
Bose -Einstein correlations. 

To conclude this section, we have shown that the sudden quench approximation in the 
evolution of the disoriented chiral condensate leads to a substantial amount of squeezing 
which manifests itself in the dramatic difference between charged and neutral pion distribu- 
tions. For an adiabatic quench the difference is much less, so that both the total multiplicity 
distributions of charged and neutral pions as well as their second order correlation functions 
are dramatic characteristic signals for the DCC and are related directly to the way in which 
the DCC forms. These are unambiguous, therefore they must be examined thoroughly in 
searches for the DCC 

6 Conclusion. 

To conclude, in this paper we have constructed effective Hamiltonians for the evolution of 
the disoriented chiral condensate without and with orientation in isospin space, starting 
from an 0(4) sigma model through the inclusion of second order quantum fluctuations. 
We have shown that both Hamiltonians have SU(1,1) symmetries leading to the presence 
of squeezed staes in their dynamics. Unlike most earlier studies, our calculations of the 
effective Hamiltonians are not restricted to zero momentum and take care of back to back 
momentum correlations. The evolution of the wave function in an expanding metric, for one 
case (without orientation in isospin space) has been considered in detail. The competing 



29 



effects of the expansion rate and the rolling down time of the system from a state of restored 
to broken symmetry has been explicitly examined. We find that in the quenched limit ( 
fast expansion) the low momentum modes are enhanced signaling DCC formation, whereas 
in the adiabatic (slow expansion) limit, no such enhancement occurs. This has been shown 
to be directly related to the value of the squeezing parameter. The manifestation of this 
difference shows up directly in the total neutral and charged multiplicity distributions at 
zero momentum and the second order correlation functions. 

Further work is required to incorporate the effects of isospin. This is being done in a 
subsequent publication [30J. Furthermore, the evolution of the wave function corresponding 
to the Hamiltonian with orientation in isospin space is also the subject of a forthcoming 
communication |31j . Within this framework back to back momentum correlations and event 
by event analysis of the experimental signals can also be done. This will provide a solid 
unified picture of the various models of DCC formation , evolution and decay, together with 
new experimental signals. Finally, we have achieved our goal of answering Bjorken's DCC 
troubling questions posed at the Trento meeting 

• Question a: Are coherent states the right quantum DCC description, or should one go 
beyond to squeezed states, etc. 

• Answer a: YES, one is naturally lead to squeezed states from a quantum field theoretic 
perspective. 

• Question b: How does one link DCC thinking to quantum effects in the data, especially 
the BoseEinstein correlations. Is DCC just another way of talking about the same 
thing? 

• Answer b: DCC formation leads to a significant difference between the charged and 
neutral pion Bose-Einstein correlations corresponding to large squeezing. So, it is not 
another way of talking about the same thing. 

• Question c: How does one enforce quantum number conservation, especially charge? 

• Answer c: Charge conservation is automatically guaranteed by the isospin analysis 
shown in paper |3U] . 

• Question d: DCC production may imply anomalous bremsstrahlung, due to the large 
quantum fluctuations in charge. Can this be calculated from first principles? 

• Answer d: This can be easily incorporated in our model and will be shown in a later 
communication. 

• Question e: Can one really set up the problem at the quantum field theory level? 

• Answer e: An emphatic YES! Albeit, the full quantum theoretical effective Hamiltonian 
with the contribution of the quarks still remains to be considered. 
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8 Appendix 1: Relationship between the transmission coefficients 
and squeezing parameter for a time dependent harmonic oscil- 
lator 

The process of particle creation and the excitation of parametric oscillators can be related 
to transmission and reflection of shrodinger waves over a potential barrier. The method we 
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describe for the relation is closely related to ref [13] Let us consider the evolution of the wave 
function of the equation 

ilk + (k 2 - V i)k (v)) ti,k = (128) 

if we define k 2 — K,&(?7) — ^{v) 2 ^ then the above equation is a time dependent harmonic 
oscillator + uj 2 (i])^(t) = . In the event of a change in frequency of the oscillator from 
co>_ to u + the asymptotic form of the real solution at 77 — ► ±00 is 

i±(t) = \{a ± e^ t + ale-^ t ) (129) 

We compare this with the complex solution of a one-d oscillator treated as a reflection 
over a barrier k 2 — with 77 corresponding to the spatial variable of a Shrodinger type of 
equation : 



£ c _(t) = e™- 1 + Re-*"-* 

£ c+ (t) = Te^' (130) 
To identify a_ = 1 + R* and a + = T 

Now consider the particle creation problem modelled by a time dependent harmonic 
oscillator [3^. The solution can be represented by 

Sr- = ^= (131) 



C r+ (0 = a — + , 
Now with the definition of the transmission coefficient 





\ 2 Vout 




2 Vifl 



T = ',° ,' oul (132) 



If we regard t as a spatial variable then these represent reflection and transmission over 
a one-d barrier. In particle creation problems the potential barrier reflection and the trans- 
mission coefficient can be related to the squeezing parameter in a particle creation problem 
by the relation Sinh 2 r = ^ =< n >= j3j3* 
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